CHAPTER (1 

) : NUMERICAL ANALYSIS 

Curve Fitting 


Method 1: Least Mean Square 
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Root Mean Square Error ( RMSE ) 
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Some Common Ideas 
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General Form 
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Interpolation : 


1- Lagrange Method 
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Newton Divided method 


P n (x) = y 0 + 8(x — x 0 ) + 8 2 (x — x„)(x — * 1 ) + 5 3 (x — x 0 )(x — x x )(x — x 2 ) 
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Operators 


1- Shift Operator E n /(x) = /(x + nh ) 

2- Forward Operator A/(x) = /(x + h ) — /(x) ; A = E — 1 

3- Backward Operator V/(x) = /(x) — /(x — h) ; V = 1 — E -1 

4- Center Operator 5 /(x) = / (x + — / (x — ; 5 — E^ — E~ 

1 _1 

5- Mean Operator nf(x) = ^ [/ (* + \) ~ / (* - fj\ = \ ^ ; 8 = ¥zl± 

A (J.g) = A/.flf r+1 +f.Ag 
A /A = gA/- Ag.f 
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Derivatives by operators form : 
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A-A 2 H—A 3 -. 
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Newton Backward 
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Gauss Forward 
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x — x 0 
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Gauss Backward 
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Stirling Form 
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Approximate for first derivatives 


a) Divided Difference Methods 


Backward divided difference 

Central divided difference 

Forward divided difference 
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Richardson Extrapolation 
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Romberg 


Therefore 
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Numerical Solution for ODE 




y' = f(x,y ) ; y(x 0 ) = y 0 ;h = — - — 

(one step) 


1- Euler Method 

yt+i =yi + hfi 

Global error 

\E n \=^[{l + hMY-l\ 

2- Taylor Method 
2 nd Order 

h 2 , 

yt+i ~ Vi + hfi + fi 

h 2 h 3 

yi+i = yi + hfi + — fi + — fi 


3 rd Order 











3- 4 th Order RK ( Runge Kutta ) 

Vi+i = y; + g ( ^i + 2 K 2 +2 K 3 + K A ) 

Ki = f( x i >Vi ) 

/ h h \ 

K i=f{* + 2 ■r + 2 K ') 

( h h \ 
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K a = f(x + h , y + h k 3 ) 

4- Modified Euler ( Corrector - Predictor) 

Predictor Step 

yf+i = y* + h fi = y£+l 

Vti i +1) = y; + ^ (fi( x i>yO + f{ x i+i 'Vi+i )) 

When you have a 2 nd Order ODE we convert it into a 1 st order system of equations then we 
solve it 
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Numerical Solution of PDE 
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Proofs 


Error by Lagrange 
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